History and current status of embryogenic culture-based tissue culture, transformation and gene editing of maize (Zea mays L.).

The production of embryogenic callus and somatic embryos is integral to the genetic improvement of crops via genetic transformation and gene editing. Regenerable embryogenic cultures also form the backbone of many micro-propagation processes for crop species. In many species, including maize, the ability to produce embryogenic cultures is highly genotype dependent. While some modern transformation and genome editing methods reduce genotype dependence, these efforts ultimately fall short of producing truly genotype-independent tissue culture methods. Recalcitrant genotypes are still identified in these genotype-flexible processes, and their presence is magnified by the stark contrast with more amenable lines, which may respond more efficiently by orders of magnitude. This review aims to describe the history of research into somatic embryogenesis, embryogenic tissue cultures, and plant transformation, with particular attention paid to maize. Contemporary research into genotype-flexible morphogenic gene-based transformation and genome engineering is also covered in this review. The rapid evolution of plant biotechnology from nascent technologies in the latter half of the 20th century to well-established, work-horse production processes has, and will continue to, fundamentally changed agriculture and plant genetics research.

A key to totipotency: Wuschel-like homeobox 2a unlocks embryogenic culture response in maize (Zea mays L.).

The ability of plant somatic cells to dedifferentiate, form somatic embryos and regenerate whole plants in vitro has been harnessed for both clonal propagation and as a key component of plant genetic engineering systems. Embryogenic culture response is significantly limited, however, by plant genotype in most species. This impedes advancements in both plant transformation-based functional genomics research and crop improvement efforts. We utilized natural variation among maize inbred lines to genetically map somatic embryo generation potential in tissue culture and identify candidate genes underlying totipotency. Using a series of maize lines derived from crosses involving the culturable parent A188 and the non-responsive parent B73, we identified a region on chromosome 3 associated with embryogenic culture response and focused on three candidate genes within the region based on genetic position and expression pattern. Two candidate genes showed no effect when ectopically expressed in B73, but the gene Wox2a was found to induce somatic embryogenesis and embryogenic callus proliferation. Transgenic B73 cells with strong constitutive expression of the B73 and A188 coding sequences of Wox2a were found to produce somatic embryos at similar frequencies, demonstrating that sufficient expression of either allele could rescue the embryogenic culture phenotype. Transgenic B73 plants were regenerated from the somatic embryos without chemical selection and no pleiotropic effects were observed in the Wox2a overexpression lines in the regenerated T0 plants or in the two independent events which produced T1 progeny. In addition to linking natural variation in tissue culture response to Wox2a, our data support the utility of Wox2a in enabling transformation of recalcitrant genotypes.

Chromosome-level genome assembly of a regenerable maize inbred line A188.

BACKGROUND: The maize inbred line A188 is an attractive model for elucidation of gene function and improvement due to its high embryogenic capacity and many contrasting traits to the first maize reference genome, B73, and other elite lines. The lack of a genome assembly of A188 limits its use as a model for functional studies. RESULTS: Here, we present a chromosome-level genome assembly of A188 using long reads and optical maps. Comparison of A188 with B73 using both whole-genome alignments and read depths from sequencing reads identify approximately 1.1 Gb of syntenic sequences as well as extensive structural variation, including a 1.8-Mb duplication containing the Gametophyte factor1 locus for unilateral cross-incompatibility, and six inversions of 0.7 Mb or greater. Increased copy number of carotenoid cleavage dioxygenase 1 (ccd1) in A188 is associated with elevated expression during seed development. High ccd1 expression in seeds together with low expression of yellow endosperm 1 (y1) reduces carotenoid accumulation, accounting for the white seed phenotype of A188. Furthermore, transcriptome and epigenome analyses reveal enhanced expression of defense pathways and altered DNA methylation patterns of the embryonic callus. CONCLUSIONS: The A188 genome assembly provides a high-resolution sequence for a complex genome species and a foundational resource for analyses of genome variation and gene function in maize. The genome, in comparison to B73, contains extensive intra-species structural variations and other genetic differences. Expression and network analyses identify discrete profiles for embryonic callus and other tissues.

Maize transformation: history, progress, and perspectives.

Maize functional genomics research and genetic improvement strategies have been greatly accelerated and refined through the development and utilization of genetic transformation systems. Maize transformation is a composite technology based on decades' efforts in optimizing multiple factors involving microbiology and physical/biochemical DNA delivery, as well as cellular and molecular biology. This review provides a historical reflection on the development of maize transformation technology including the early failures and successful milestones. It also provides a current perspective on the understanding of tissue culture responses and their impact on plant regeneration, the pros and cons of different DNA delivery methods, the identification of a palette of selectable/screenable markers, and most recently the development of growth-stimulating or morphogenic genes to improve efficiencies and extend the range of transformable genotypes. Steady research progress in these interdependent components has been punctuated by benchmark reports celebrating the progress in maize transformation, which invariably relied on a large volume of supporting research that contributed to each step and to the current state of the art. The recent explosive use of CRISPR/Cas9-mediated genome editing has heightened the demand for higher transformation efficiencies, especially for important inbreds, to support increasingly sophisticated and complicated genomic modifications, in a manner that is widely accessible. These trends place an urgent demand on taking maize transformation to the next level, presaging a new generation of improvements on the horizon. Once realized, we anticipate a near-future where readily accessible, genotype-independent maize transformation, together with advanced genomics, genome editing, and accelerated breeding, will contribute to world agriculture and global food security.

Assessing the Viability of Recovery of Hydroxycinnamic Acids from Lignocellulosic Biorefinery Alkaline Pretreatment Waste Streams.

Invited for this month's cover is the research team from the D.O.E. Great Lake Bioenergy Research Center (GLBRC) at the University of Wisconsin-Madison. The cover image shows how a diverse team with expertise in many different fields works together in an integrated fashion to address complex problems. Only when the whole system, from field to the liquid fuels and co-products, is assessed, can we identify the key parameters needed to design an economically viable biorefinery-based economy. Cover art by Chelsea Mamott. The Full Paper itself is available at 10.1002/cssc.201903345.

Assessing the Viability of Recovery of Hydroxycinnamic Acids from Lignocellulosic Biorefinery Alkaline Pretreatment Waste Streams.

The hydroxycinnamic acids p-coumaric acid (pCA) and ferulic acid (FA) add diversity to the portfolio of products produced by using grass-fed lignocellulosic biorefineries. The level of lignin-bound pCA in Zea mays was modified by the alteration of p-coumaroyl-CoA monolignol transferase expression. The biomass was processed in a lab-scale alkaline-pretreatment biorefinery process and the data were used for a baseline technoeconomic analysis to determine where to direct future research efforts to couple plant design to biomass utilization processes. It is concluded that future plant engineering efforts should focus on strategies that ramp up accumulation of one type of hydroxycinnamate (pCA or FA) predominantly and suppress that of the other. Technoeconomic analysis indicates that target extraction titers of one hydroxycinnamic acid need to be >50 g kg-1 biomass, at least five times higher than observed titers for the impure pCA/FA product mixture from wild-type maize. The technical challenge for process engineers is to develop a viable process that requires more than 80 % reduction of the isolation costs.

Maize sugary enhancer1 (se1) is a gene affecting endosperm starch metabolism.

sugary enhancer1 (se1) is a naturally occurring mutant allele involved in starch metabolism in maize endosperm. It is a recessive modifier of sugary1 (su1) and commercially important in modern sweet corn breeding, but its molecular identity and mode of action remain unknown. Here, we developed a pair of near-isogenic lines, W822Gse (su1-ref/su1-ref se1/se1) and W822GSe (su1-ref/su1-ref Se1/Se1), that Mendelize the se1 phenotype in an su1-ref background. W822Gse kernels have lower starch and higher water soluble polysaccharide and sugars than W822GSe kernels. Using high-resolution genetic mapping, we found that wild-type Se1 is a gene Zm00001d007657 on chromosome 2 and a deletion of this gene causes the se1 phenotype. Comparative metabolic profiling of seed tissue between these 2 isolines revealed the remarkable difference in carbohydrate metabolism, with sucrose and maltose highly accumulated in the mutant. Se1 is predominantly expressed in the endosperm, with low expression in leaf and root tissues. Differential expression analysis identified genes enriched in both starch biosynthesis and degradation processes, indicating a pleiotropic regulatory effect of se1 Repressed expression of Se1 and Su1 in RNA interference-mediated transgenic maize validates that deletion of the gene identified as Se1 is a true causal gene responsible for the se1 phenotype. The findings contribute to our understanding of starch metabolism in cereal crops.

Genome-wide association analysis of stalk biomass and anatomical traits in maize.

BACKGROUND: Maize stover is an important source of crop residues and a promising sustainable energy source in the United States. Stalk is the main component of stover, representing about half of stover dry weight. Characterization of genetic determinants of stalk traits provide a foundation to optimize maize stover as a biofuel feedstock. We investigated maize natural genetic variation in genome-wide association studies (GWAS) to detect candidate genes associated with traits related to stalk biomass (stalk diameter and plant height) and stalk anatomy (rind thickness, vascular bundle density and area). RESULTS: Using a panel of 942 diverse inbred lines, 899,784 RNA-Seq derived single nucleotide polymorphism (SNP) markers were identified. Stalk traits were measured on 800 members of the panel in replicated field trials across years. GWAS revealed 16 candidate genes associated with four stalk traits. Most of the detected candidate genes were involved in fundamental cellular functions, such as regulation of gene expression and cell cycle progression. Two of the regulatory genes (Zmm22 and an ortholog of Fpa) that were associated with plant height were previously shown to be involved in regulating the vegetative to floral transition. The association of Zmm22 with plant height was confirmed using a transgenic approach. Transgenic lines with increased expression of Zmm22 showed a significant decrease in plant height as well as tassel branch number, indicating a pleiotropic effect of Zmm22. CONCLUSION: Substantial heritable variation was observed in the association panel for stalk traits, indicating a large potential for improving useful stalk traits in breeding programs. Genome-wide association analyses detected several candidate genes associated with multiple traits, suggesting common regulatory elements underlie various stalk traits. Results of this study provide insights into the genetic control of maize stalk anatomy and biomass.

Genetic Fine-Mapping of a Quantitative Trait Locus (QTL) Associated with Embryogenic Tissue Culture Response and Plant Regeneration Ability in Maize (Zea mays L.).

Embryogenic and regenerable tissue cultures are widely utilized in plant transformation, clonal propagation, and biological research applications. Germplasm utilized in those applications are limited, however, due to genotype-dependent culture response. The goal of this study was to identify genomic regions controlling embryogenic and regenerable tissue culture response in the globally important crop, maize ( L.), toward the long-term objective of developing approaches for genotype-independent plant genetic engineering and clonal propagation systems. An inbred maize line, WCIC2, nearly-isogenic to reference inbred B73, was developed by phenotypic selection and molecular marker analysis. WCIC2 has over 50x increase in tissue culture response relative to the recurrent parent, B73. This line was used to genetically fine-map a region on chromosome 3 controlling embryogenic and regenerable tissue culture response to a 23.9 Mb region. WCIC2 and derivatives will be useful materials to enable maize research in a genetic background similar to B73, and our genetic mapping results will advance research to identify causal genes controlling somatic embryo formation and plant regeneration in maize.

Suppression of CINNAMOYL-CoA REDUCTASE increases the level of monolignol ferulates incorporated into maize lignins.

BACKGROUND: The cell wall polymer lignin provides structural support and rigidity to plant cell walls, and therefore to the plant body. However, the recalcitrance associated with lignin impedes the extraction of polysaccharides from the cell wall to make plant-based biofuels and biomaterials. The cell wall digestibility can be improved by introducing labile ester bonds into the lignin backbone that can be easily broken under mild base treatment at room temperature. The FERULOYL-CoA MONOLIGNOL TRANSFERASE (FMT) enzyme, which may be naturally found in many plants, uses feruloyl-CoA and monolignols to synthesize the ester-linked monolignol ferulate conjugates. A mutation in the first lignin-specific biosynthetic enzyme, CINNAMOYL-CoA REDUCTASE (CCR), results in an increase in the intracellular pool of feruloyl-CoA. RESULTS: Maize (Zea mays) has a native putative FMT enzyme, and its ccr mutants produce an increased pool of feruloyl-CoA that can be used for conversion to monolignol ferulate conjugates. The decreased lignin content and monomers did not, however, impact the plant growth or biomass. The increase in monolignol conjugates correlated with an improvement in the digestibility of maize stem rind tissue. CONCLUSIONS: Together, increased monolignol ferulates and improved digestibility in ccr1 mutant plants suggests that they may be superior biofuel crops.

Autophagic recycling plays a central role in maize nitrogen remobilization.

Autophagy is a primary route for nutrient recycling in plants by which superfluous or damaged cytoplasmic material and organelles are encapsulated and delivered to the vacuole for breakdown. Central to autophagy is a conjugation pathway that attaches AUTOPHAGY-RELATED8 (ATG8) to phosphatidylethanolamine, which then coats emerging autophagic membranes and helps with cargo recruitment, vesicle enclosure, and subsequent vesicle docking with the tonoplast. A key component in ATG8 function is ATG12, which promotes lipidation upon its attachment to ATG5. Here, we fully defined the maize (Zea mays) ATG system transcriptionally and characterized it genetically through atg12 mutants that block ATG8 modification. atg12 plants have compromised autophagic transport as determined by localization of a YFP-ATG8 reporter and its vacuolar cleavage during nitrogen or fixed-carbon starvation. Phenotypic analyses showed that atg12 plants are phenotypically normal and fertile when grown under nutrient-rich conditions. However, when nitrogen-starved, seedling growth is severely arrested, and as the plants mature, they show enhanced leaf senescence and stunted ear development. Nitrogen partitioning studies revealed that remobilization is impaired in atg12 plants, which significantly decreases seed yield and nitrogen-harvest index. Together, our studies demonstrate that autophagy, while nonessential, becomes critical during nitrogen stress and severely impacts maize productivity under suboptimal field conditions.

Whole transcriptome profiling of maize during early somatic embryogenesis reveals altered expression of stress factors and embryogenesis-related genes.

Embryogenic tissue culture systems are utilized in propagation and genetic engineering of crop plants, but applications are limited by genotype-dependent culture response. To date, few genes necessary for embryogenic callus formation have been identified or characterized. The goal of this research was to enhance our understanding of gene expression during maize embryogenic tissue culture initiation. In this study, we highlight the expression of candidate genes that have been previously regarded in the literature as having important roles in somatic embryogenesis. We utilized RNA based sequencing (RNA-seq) to characterize the transcriptome of immature embryo explants of the highly embryogenic and regenerable maize genotype A188 at 0, 24, 36, 48, and 72 hours after placement of explants on tissue culture initiation medium. Genes annotated as functioning in stress response, such as glutathione-S-transferases and germin-like proteins, and genes involved with hormone transport, such as PINFORMED, increased in expression over 8-fold in the study. Maize genes with high sequence similarity to genes previously described in the initiation of embryogenic cultures, such as transcription factors BABY BOOM, LEAFY COTYLEDON, and AGAMOUS, and important receptor-like kinases such as SOMATIC EMBRYOGENESIS RECEPTOR LIKE KINASES and CLAVATA, were also expressed in this time course study. By combining results from whole genome transcriptome analysis with an in depth review of key genes that play a role in the onset of embryogenesis, we propose a model of coordinated expression of somatic embryogenesis-related genes, providing an improved understanding of genomic factors involved in the early steps of embryogenic culture initiation in maize and other plant species.

Interaction of photoperiod and vernalization determines flowering time of Brachypodium distachyon.

Timing of flowering is key to the reproductive success of many plants. In temperate climates, flowering is often coordinated with seasonal environmental cues such as temperature and photoperiod. Vernalization is an example of temperature influencing the timing of flowering and is defined as the process by which a prolonged exposure to the cold of winter results in competence to flower during the following spring. In cereals, three genes (VERNALIZATION1 [VRN1], VRN2, and FLOWERING LOCUS T [FT]) have been identified that influence the vernalization requirement and are thought to form a regulatory loop to control the timing of flowering. Here, we characterize natural variation in the vernalization and photoperiod responses in Brachypodium distachyon, a small temperate grass related to wheat (Triticum aestivum) and barley (Hordeum vulgare). Brachypodium spp. accessions display a wide range of flowering responses to different photoperiods and lengths of vernalization. In addition, we characterize the expression patterns of the closest homologs of VRN1, VRN2 (VRN2-like [BdVRN2L]), and FT before, during, and after cold exposure as well as in different photoperiods. FT messenger RNA levels generally correlate with flowering time among accessions grown in different photoperiods, and FT is more highly expressed in vernalized plants after cold. VRN1 is induced by cold in leaves and remains high following vernalization. Plants overexpressing VRN1 or FT flower rapidly in the absence of vernalization, and plants overexpressing VRN1 exhibit lower BdVRN2L levels. Interestingly, BdVRN2L is induced during cold, which is a difference in the behavior of BdVRN2L compared with wheat VRN2 during cold.

Assessing the efficiency of RNA interference for maize functional genomics.

A large-scale functional genomics project was initiated to study the function of chromatin-related genes in maize (Zea mays). Transgenic lines containing short gene segments in inverted repeat orientation designed to reduce expression of target genes by RNA interference (RNAi) were isolated, propagated, and analyzed in a variety of assays. Analysis of the selectable marker expression over multiple generations revealed that most transgenes were transmitted faithfully, whereas some displayed reduced transmission or transgene silencing. A range of target-gene silencing efficiencies, from nondetectable silencing to nearly complete silencing, was revealed by semiquantitative reverse transcription-PCR analysis of transcript abundance for the target gene. In some cases, the RNAi construct was able to cause a reduction in the steady-state RNA levels of not only the target gene, but also another closely related gene. Correlation of silencing efficiency with expression level of the target gene and sequence features of the inverted repeat did not reveal any factors capable of predicting the silencing success of a particular RNAi-inducing construct. The frequencies of success of this large-scale project in maize, together with parameters for optimization at various steps, should serve as a useful framework for designing future RNAi-based functional genomics projects in crop plants.

The Lem2 gene promoter of barley directs cell- and development-specific expression of gfp in transgenic plants.

Transgenic approaches to combating fungal pathogens, such as Fusarium graminearum, require the targeting of antifungal gene expression in tissues of developing seed spikes of cereal grains, especially lemmas and epicarps. The Lem2 gene of barley encodes a lectin-like protein that is strongly up-regulated by salicylic acid and is preferentially expressed in lemmas, paleas (lemma/palea) and coleoptiles. Transient expression studies have indicated that the proximal -75/+70 region (relative to the transcription start site) determines organ specificity. In the present study, Golden Promise barley stably transformed with Morex Lem2 promoter/gfp reporter constructs displayed cell- and development-specific expression of gfp (green fluorescent protein gene). This expression corresponded to the expression seen in Northern blots of Morex organs. Under the full-length promoter, strong GFP fluorescence was observed in the lemma/palea, glumes, coleoptile, auricle and ligule. Weak GFP fluorescence was also observed in the rachis, tips of primary leaves and the leaf sheath. Unexpectedly, strong expression occurred in the epicarp, even though Lem2 is not expressed in this organ in Morex. Studies showed that the Lem2 promoter is more highly methylated in the epicarp than in the lemma of Morex. In the lemma/palea, gfp underwent a temporal shift in expression from the mesophyll to specialized epidermal cork cells. Similar to the lemma/palea, expression in the leaf sheath was localized in the cork cells. Progressive 5' deletions of the promoter to nucleotide -75 gradually reduced the level of gfp expression, but tissue- and cell-specific expression was retained.

Spatial and temporal divergence of expression in duplicated barley germin-like protein-encoding genes.

Subfunctionalization is the process by which a pair of duplicated genes, or paralogs, experiences a reduction of individual expression patterns or function while still reproducing the complete expression pattern and function of the ancestral gene. Two germin-like protein (GLP)-encoding genes, GerB and GerF, are paralogs that belong to a small gene family in barley (Hordeum vulgare). Both genes share high nucleotide sequence similarity in coding and noncoding regions and encode identical apoplastic proteins. The use of RNA gel blots, coupled with single-stranded conformation polymorphism (SSCP) analysis of RT-PCR products, elucidated the developmental and tissue-specific expression patterns of each gene. Individual expression patterns provided evidence of both overlapping redundancy and early subfunctionalization. GerB is predominantly expressed in developing shoots, while GerF is predominantly expressed in seedling roots, developing spikes, and pericarp/testa. GerF promoter deletion studies located a region (-356/-97) responsible for high promoter activity and showed the ability of GerB and GerF upstream regions to drive gfp expression in coleoptiles, epicarps, and lemma/palea of developing spikes. The observed expression patterns are consistent with proposed roles in plant development and defense mechanisms for this gene family. These roles may explain why redundancy has been selectively maintained in this duplicate gene pair.

The complex developmental expression of a novel stress-responsive barley Ltp gene is determined by a shortened promoter sequence.

The search for a cereal promoter capable of driving preferential transgene expression in the pericarp epidermis (epicarp) of developing barley (Hordeum vulgare L.) resulted in the cloning of a novel gene. This encoded a polypeptide of 124 amino acids showing 87 identity with WBP1A, a wheat lipid transfer protein (LTP), but much lower homology to other barley LTPs. In addition to the epicarp, this Ltp-like gene, Ltp6, is highly expressed in coleoptiles and embryos under normal growth conditions. Messenger RNA levels increased in seedling tissues during salt and cold treatments and under applied abscisic acid (ABA) and salicylic acid (SA). Taken together, Ltp6 tissue-specific and response patterns are distinct from other known barley Ltp genes. Inverse PCR was used to derive 2345 bp of upstream Ltp6 sequence. The level of transcription conferred by different promoter deletion constructs was assessed by quantitative real time RT-PCR using gfp as a reporter in transient expression assays. All constructs containing at least 192 bp of upstream sequence and the 5'UTR conferred tissue-specific expression and retained most of the promoter strength. Deletion of 64 bp (-192/-128) from this upstream sequence reduced expression levels by 80. Moreover, a minimal 247 bp Ltp6 promoter continuously drove gfp expression during spike development, from early ovary differentiation through its final expression in the epicarp and during embryogenesis and germination in transgenic barley, reproducing the expression pattern of the native gene. The potential use of this promoter sequence for targeting transgene-mediated disease resistance in barley and wheat is discussed.

Transgene-induced RNA interference as a tool for plant functional genomics.

RNA interference (RNAi) is a powerful tool for functional genomics in a number of species. The logistics and procedures for doing high-throughput RNAi to investigate the functions of large numbers of genes in Arabidopsis thaliana and in Zea mays are described. Publicly available plasmid vectors that facilitate the stable chromosomal integration of inverted repeat transgenes that trigger RNAi have been used to generate more than 50 independent transgenic lines each in Arabidopsis and maize. Analysis of mRNA abundance of the targeted genes in independent lines transformed with distinct constructs indicates that the success of RNAi-induced silencing is gene dependent. mRNA levels were not detectably reduced for some genes, but were dramatically reduced for a number of genes targeted. A common pattern was that multiple independent lines transgenic for the same construct showed the same extent of silencing. This chapter describes the procedures used to generate and test transgenic lines mediating RNAi in Arabidopsis and maize.

A proximal upstream sequence controls tissue-specific expression of Lem2, a salicylate-inducible barley lectin-like gene.

The lemma and palea (lemma/palea), which form the husk of barley (Hordeum vulgare L.) seeds, constitutively express high levels of defense-related genes, relative to leaves [Abebe et al. (2004) Crop Sci 44:942-950]. One of these genes, Lem2, is expressed mainly in the lemma/palea and coleoptile and is strongly upregulated by salicylic acid (SA) and its functional analog 2,6-dichloroisonicotinic acid . Induction by SA was rapid, occurring within 4 h of treatment. However, Lem2 is not responsive to methyl jasmonate (MeJA) or wounding and is downregulated by drought, dehydration, and abscisic acid. These results suggest that Lem2 is involved in systemic acquired resistance. Sequence analysis showed that LEM2 is a jacalin-related lectin (JRL)-like protein with two domains. Consistent with northern and western blot data, transient expression analyses using Lem2::gfp constructs showed strong expression in lemmas and a trace expression in leaves. Successive 5' deletions of the 1,414 bp upstream region gradually weakened promoter strength, as measured by real-time PCR. Promoter deletion studies also revealed that the -75/+70 region (containing the TATA box, 5' UTR, and a SA-response element) determines tissue specificity and that the distal promoter region simply enhances expression. Southern analysis indicated that Morex barley has at least three copies of the Lem2 gene arranged in tandem on chromosome 5(1H) Bin 02, near the short arm telomere. Lem2 is not present in the barley cultivars Steptoe, Harrington, Golden Promise, and Q21861.

Comparative analysis of SET domain proteins in maize and Arabidopsis reveals multiple duplications preceding the divergence of monocots and dicots.

Histone proteins play a central role in chromatin packaging, and modification of histones is associated with chromatin accessibility. SET domain [Su(var)3-9, Enhancer-of-zeste, Trithorax] proteins are one class of proteins that have been implicated in regulating gene expression through histone methylation. The relationships of 22 SET domain proteins from maize (Zea mays) and 32 SET domain proteins from Arabidopsis were evaluated by phylogenetic analysis and domain organization. Our analysis reveals five classes of SET domain proteins in plants that can be further divided into 19 orthology groups. In some cases, such as the Enhancer of zeste-like and trithorax-like proteins, plants and animals contain homologous proteins with a similar organization of domains outside of the SET domain. However, a majority of plant SET domain proteins do not have an animal homolog with similar domain organization, suggesting that plants have unique mechanisms to establish and maintain chromatin states. Although the domains present in plant and animal SET domain proteins often differ, the domains found in the plant proteins have been generally implicated in protein-protein interactions, indicating that most SET domain proteins operate in complexes. Combined analysis of the maize and Arabidopsis SET domain proteins reveals that duplication of SET domain proteins in plants is extensive and has occurred via multiple mechanisms that preceded the divergence of monocots and dicots.